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Abstract 

We applied a mean field approach to study clustering coefficients in Barabasi- Albert networks. 
We found that the local clustering in BA networks depends on the node degree. Analytic results 
have been compared to extensive numerical simulations finding a very good agreement for nodes 
with low degrees. Clustering coefficient of a whole network calculated from our approach perfectly 
fits numerical data. 

PACS numbers: 89.75.-k, 02.50.-r, 05.50.+q 
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Introduction. During the last decade networks became a very popular research domain 
among physicists (for a review see . It is not surprising, since networks are every- 

where. They surround us. In our daily life we participate in dozens of them. A number of 
social institutions, communication and biological systems may be represented as networks 
i.e. sets of nodes connected by links. It was observed that despite functional diversity most 
of real web-like systems share similar structural properties. The properties are: fat-tailed 
degree distribution (that allows for hubs i.e. nodes with high degree), small average distance 
between any two nodes (the so-called small world effect) and a large penchant for creating 
chques (i.e. highly interconnected groups of nodes). 

A number of network construction procedures have been proposed to incorporate the 
characteristics. The Barabasi- Albert (BA) [J, |^ growing network model is probably the 
best known. Two important ingredients of the model are: continuous network growth and 
preferential attachment. The network starts to grow from an initial cluster of m fully 
connected sites. Each new node that is added to the network creates m links that connect 
it to previously added nodes. The preferential attachment means that the probability of a 
new link to end up in a vertex i is proportional to the connectivity ki of this vertex 

n, = m-^. (1) 

Taking into account that J2jkj = 2mt the last formula may be rewritten as Ilj = fcj/(2t). 
By means of mean field approximation one can find that the average degree of a node i 
that entered the network at time ti increases with time as a power-law 

h{t)=mJj. (2) 
V ti 

Taking advantage of the above formula one can calculate the probability that two randomly 
selected nodes i and j are nearest neighbors. It is given by 

m 1 

It was shown that the degree distribution in BA network follows a power-law 

m = 9-. (4) 

where A; = m, m -|- 1, . . . , m\/t. The power law degree distribution is characteristic of many 
real- world networks and the scaling exponent as a = 3 is close to those observed in real 



systems areai £ (2,3). It was also shown that the BA model is a small world. The mean 
distance between sites in the network having t nodes behaves as I ~ In t/ In In t 0, Q|- The 
only shortcoming of the model is that it does not incorporate a high degree of cliqueness 
observed in real networks. 

In this paper we study cliqueness effects in BA networks. The cliqueness is measured by 
the clustering coefficient C jsj, 0]. The clustering coefficient Ci of a single node i describes 
the density of connections in the neighborhood of this node. It is given by the ratio of the 
number Ei of links between the nearest neighbors of i and the potential number of such links 

Emax = ki{ki l)/2 

'~ E ~ k(k- - IV ^ ' 

-'-'max -'-J 

The clustering coefficient C of the whole network is the average of all individual Cj's. It is 
known, from numerical calculations, that the clustering coefficient in BA networks rapidly 
decreases with the network size t. In this article we apply a mean field approach to study the 
parameter. Our calculations confirm that in the limit of large {t ^ 1) and dense (m ^ 1) 
. 3A networks the clustering coefficient scales as the clustering coefficient in random graphs 
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with an appropriate scale-free degree distribution 



8 t ^ ^ 



We also show that the individual clustering coefficient Ci in BA network weakly depends 
on node's degree fcj. The dependence is almost invisible when one looks at numerical data 
presented by other authors ,13]. 

Mean field approach. Let us concentrate on a certain node i in a BA network of size t. 
We assume that m > 2. The case of m = 1 is trivial. BA networks with m = 1 are trees 
thus the clustering coefficient in these networks is equal to zero. By the definition (0) the 
clustering coefficient Ci depends on two variables Ei and fcj. Since in the BA model only 
new nodes may create links the coefficient Ci changes only when its degree ki changes i.e. 
when new nodes create connections to i and x G (0,m — 1) of its nearest neighbors. The 
appropriate equation for changes of Ci is then 

m-l 

-f^ =Y1 Pix^Ci^ , (7) 

where /S.Cix denotes the change of the clustering coefficient when a new node connects to the 
node i and to x of the first neighbors of z, whereas pix describes the probability of this event. 



ACix is simply the difference between clustering coefficients of the same node i calculated 
after and before a new node attachment 

2{E, + x) 2E, 2C, 2x 

" h{h + 1) Uh - 1) A;, + 1 ^ UK + 1) ■ ^ ' 

The probability is a product of two factors. The first factor is the probability of a new 

link to end up in i. The probability is given by (0). The second one is the probability 

that among the rest of (m — 1) new links x links connect to nearest neighbors of i. It is 

equivalent to the probability that (m — 1) Bernoulli trials with the probability for success 

equal to Y,j*kj / Hv^v = J2j*kj/{2mt) result in x successes runs over the nearest 

neighbors of the node i). Replacing the sum by an integral one obtains 

/t 
kj Pij dtj = — kilnt. (9) 
J* 

Summarizing the above discussion one yields the relation 

Now, inserting ((HI) and (fTUj) into ((7j) one obtains after some algebra 

dCi m m(m — 1) In t , . 

* + 77 — ,9 , rrr\- \^^) 



dt {mt + A{mt'^ + 

Solving the equation for Ci one gets 

cm = ^'^-}} { {Intr - l./^lnt - -^I^ + b] , (12) 

where B is an integration constant and may be determined from the initial condition Ci{ti) 
that describes the clustering coefficient of the node i exactly at the moment of its attachment 
U 

aw4Ei:P«P»ft./(")-5^^. (13) 

Following the notation introduced by Bianconi and Capocci the initial clustering coef- 
ficient Ciiti) may may be written as 



Ci{ti) — 



dt 



(14) 

t=ti 

where d{N'i{t)) / dt describes how the number of triangular loops increases in time. FigC] 
shows the prediction of the equation ()13|) in comparison with numerical results. For small 
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FIG. 1: The initial value of the local clustering coefficient Cj(tj) (averaged over 1000 BA networks). 

values of ti the numerical data differ from the theory in a significant way. This can be 
explained by the fact that the formula for the probability of a connection pij (jS)), that we 



Taking into account the initial condition Ci{ti) and neglecting mutually compensating 
terms that occure in (fT^ after putting B calculated from ()13|) one obtains the formula for 
time evolution of the clustering coefficient of a given node i 



Let us note that if ^ t or m ^ 1 the local clustering coefficient does not depend on the 
node under consideration and approaches Ci{t) ~ (m — l)(lnt)^/(8t) i.e. the formula © that 
gives the the clustering coefficient of a random graph with a power-low degree distribution 
(0}. Since one knows how the node's degree evolves in time one can also calculate the 
formula for Ci{k.j). At the Fig|21we present the formula (fT^ (solid line) and corresponding 
numerical data (scatter plots). The two kinds of scatter plots represent respectively: real 
data (light gray circles) and the same data subjected adjacent averaging smoothing (dark 
gray circles). As before (see Fig.([Q)), we observe a significant disagreement between the 
numerical data and the theory for small ti. We suspect that the effect has the same origin 
i.e. the relations (j2I) and (jH)) that we use in our derivation work well only in the asymptotic 
region <t — > oo. 

To obtain the clustering coefficient C of the whole network the expression (fT^ has to 



use three times in (jl3|) . holds only in the asymptotic region tj 



oo. 




(15) 
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FIG. 2: The local clustering coefficient Cj(t) as a function tj (averaged over 10^ networks). Note 
that the ki axis is nonlinear. 

be averaged over all nodes within a network C = Jl Ci{t)dti/t. We were not able to find 
an exact analytic form of this integral but corresponding numerical values (open squares at 
the Fig|2)) fit very well a mean field approximation that we propose below (solid line at the 
FigOl) 

f^_l)((lnt)2 + ^^(lnt,)2\ 

' 8 (iVt + Vtl/my)^ ^ ^ 

After performing separate integration of the numerator and the denominator one gets 

^ _ ((m + l)2(lnt)2 - 8mlnt + 8m) 

~ 8(m- l)(6m2 + 8m + 3)t ' ^ ^ 

For large (t — > oo) and dense (m ^ 1) networks the above formula approaches ©. 
The effect lets us deduce that the structural correlations T5] characteristic for growing BA 
networks become less important in larger and denser networks. The same was suggested 
in . FigO shows the average clustering coefficient in BA networks as a function of the 
network size t compared with the analytical formula (jl7p . 

Conclusions. In summary, we applied a mean field approach to study clustering effects in 
Barabasi- Albert networks. We found that the BA networks do not show the homogeneous 
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lisj ]. We derived a general formula for the clustering coeffi- 
cient C characterizing the whole BA network. We found that in the limit of large [t — > oo) 



clustering as suggested in 
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FIG. 3: The clustering coefficient C of a whole BA network as a function of the network size t 
(averaged over 100 networks). 

and dense (m ^ 1) networks both the local (Cj) and the global (C) clustering coefficients 
approach clustering coefficient derived for a random graph with a power-low degree distribu- 
tion (jH). Our derivations were checked against numerical simulation of BA networks finding 
a very good agreement. 
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